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We study the classical dynamics of resonantly modulated large-spin systems in a strong magnetic 
field. We show that these systems have special symmetry. It leads to characteristic nonlinear 
effects. They include abrupt switching between magnetization branches with varying modulating 
field without hysteresis and a specific pattern of switching in the presence of multistability and 
hysteresis. Along with steady forced vibrations the transverse spin components can display transient 
vibrations at a combination of the Larmor frequency and a slower frequency determined by the 
anisotropy energy. The analysis is based on a microscopic theory that takes into account relaxation 
mechanisms important for single-molecule magnets and other large-spin systems. We find how the 
Landau-Lifshitz model should be modified in order to describe the classical spin dynamics. The 
occurrence of transient oscillations depends on the interrelation between the relaxation parameters. 
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I. INTRODUCTION 

Large-spin systems have a finite but comparatively 
large number of quantum states. Therefore a single sys- 
tem can be used to study a broad range of phenomena, 
from purely quantum to semiclassical where the spin be- 
haves almost like a classical top. One of the interesting 
features of large-spin systems is that, in a strong static 
magnetic field, their energy levels become almost equidis- 
tant, with level spacing close to hu>o, where loq is the Lar- 
mor frequency. As a result, radiation at frequency ~ u>o is 
resonant simultaneously for many inter level transitions. 
This leads to new quantum and classical nonlinear reso- 
nant effects. 

An important class of large-spin systems is single- 
molecule magnets (SMMs). SMMs display an extremely 
rich behavior and have been attracting much attention 
in recent years. A variety of SMMs has already been 
discovered and investigated, see Refs. HUM! for a review, 
and new systems are being found 4 ' 5 . Another example of 
large-spin systems is provided by large nuclear spins, the 
interest in which has renewed in view of their possible 
use in quantum computing 6 . 

In this paper we study the dynamics of large-spin sys- 
tems, 5 3> 1, in the classical limit. We assume that the 
system is in a strong static magnetic field along the easy 
magnetization axis and in an almost resonant transverse 
field. For a small relaxation rate, even a weak transverse 
resonant field can lead to hysteresis of the response. As 
we show, the hysteresis is quite unusual. 

It is convenient to describe the dynamics of a reso- 
nantly modulated spin in the rotating wave approxima- 
tion (RWA). The corresponding analysis in the absence 
of relaxation has revealed a special quantum feature, an 
antiresonance of the response which accompanies anti- 
crossing of quasienergy levels^. Quantum spin dynamics 
in the rotating frame bears also on the dynamics of the 
Lipkin-Meshkov-Glick model2i2.I2.ii. 

One may expect that the features of the coherent quan- 



tum dynamics should have counterparts in the classical 
spin dynamics in the presence of dissipation. As we show, 
this is indeed the case. The system displays an unusual 
behavior in a certain range of modulation parameters. 
This behavior is due to a special symmetry. It leads to 
specific features of hysteresis and to discontinuous (in 
the neglect of fluctuations) switching between different 
response branches even in the absence of hysteresis. 

Classical dynamics of a large-spin system in a resonant 
field would be expected to have similarities with the dy- 
namics of a modulated magnetic nanoparticle near fer- 
romagnetic resonance. It was understood back in the 
1950' s 12 ' 13 that the response near ferromagnetic reso- 
nance becomes strongly nonlinear already for compar- 
atively weak radiation strength due to the magnetization 
dependence of the effective magnetic field. Resonant re- 
sponse may become multivalued as a function of the mod- 
ulating field amplitud o 14 ' 15 . A detailed analysis of non- 
linear magnetization dynamics in uniaxial nanoparticles 
modulated by a strong circularly polarized periodic field 
was done recently^. These studies as well as many other 
studies of magnetization dynamics in ferromagnets were 
based on the Landau-Lifshitz-Gilbert equation. 

In contrast to magnetic nanoparticles, for large-spin 
systems quantum effects are substantial. A distinction 
which remains important in the classical limit concerns 
relaxation mechanisms. Spin relaxation occurs via tran- 
sitions between discrete energy levels with emission, ab- 
sorption, or inelastic scattering of excitations of a thermal 
reservoir to which the spin is coupled. Relevant relax- 
ation mechanisms depend on the specific system but as 
we show, even in the classical limit relaxation is not de- 
scribed, generally, by the Landau-Lifshitz damping. As 
a result the classical spin dynamics strongly differs from 
the dynamics of a magnetic nanoparticle. 

The microscopic analysis of relaxation is simplified in 
the presence of a strong static magnetic field. Here, all 
spin energy levels are almost equidistant. Therefore exci- 
tations of the thermal bath emitted, for example, in tran- 
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sitions within different pairs of neighboring levels have 
almost the same energies. As a consequence, relaxation 
is described by a small number of constants independent 
of the form of the weighted with the interaction density 
of states of the bath, and the analysis applies for an ar- 
bitrary ratio between the level nonequidistance and their 
relaxational broadening 17 . 

We consider three relaxation mechanisms. Two of 
them correspond to transitions between neighboring and 
next neighboring spin levels, with the coupling to bosonic 
excitations quadratic in the spin operators. Such cou- 
pling is important, in particular, for SMMs where energy 
relaxation is due to phonon scattering. The theory of 
relaxation of SMMs was developed earlie r 18 ' 19 and has 
been tested experimentally, see Refs. I20H21I and papers 
cited therein. We also consider coupling to a bosonic 
bath linear in spin operators. It leads to relaxation that 
in the classical limit has the form of the Landau-Lifshitz 
damping provided the modulation field is weak compared 
to the static field, as assumed in the RWA. 

The typical duration of scattering events that lead to 
spin relaxation is often ~ . ^ n ^ ne RWA they appear 
instantaneous. The operator that describes relaxation 
has a simple functional form, with no retardation in the 
"slow" time. This is advantageous for studying the clas- 
sical limit and allows us to obtain analytical results. 

In the classical limit, a spin is characterized by two 
dynamical variables, for example, azimuthal and polar 
angles. In the RWA, they satisfy autonomous equations 
of motion, the coefficients in these equations do not de- 
pend on time. A two-variable nonlinear dissipative sys- 
tem can have both stationary and periodic states^ 2 -. As 
we show, such states indeed emerge for a resonantly mod- 
ulated spin. They were predicted also for a magnetic 
nanoparticle with Landau-Lifshitz damping in the case 
of a generally nonresonant modulation^. 

For a spin, the occurrence of periodic states in the ro- 
tating frame critically depends on the interrelation be- 
tween the relaxation parameters. In particular, these 
states do not emerge for a resonantly modulated spin 
with microscopic relaxation that reduces to the Landau- 
Lifshitz damping in the RWA. Moreover, quantum fluc- 
tuations lead to phase diffusion which results in decay of 
periodic states in the rotating frame, making the corre- 
sponding vibrations transient. 

The paper is organized as follows. In Sec. [TT| we in- 
troduce a model of the spin and its interaction with a 
thermal bath and derive the quantum kinetic equation 
with account taken of different relaxation mechanisms. 
In Sec. IIIII we obtain classical equations of motion and 
discuss the symmetry of the system. We find analyti- 
cally, for weak damping, the positions of the bifurcation 
curves where the number of stationary states in the rotat- 
ing frame changes (saddle-node bifurcations) and where 
periodic states are split off from stationary states (Hopf 
bifurcations). Sec. |IV] describes the specific and, per- 
haps, most unusual feature of the system, the occurrence 
of Hamiltonian-like dynamics in the presence of dissipa- 



tion. In Sec. E] spin dynamics and hysteresis are de- 
scribed for the relation between relaxation parameters 
where the system does not have periodic states in the 
rotating frame. In Sec. I VII we consider the opposite case. 
The onset of periodic states and their stability are an- 
alyzed and the features of the hysteresis related to the 
occurrence of periodic states are studied. Details of the 
calculations are outlined in Appendix. Sec. IVIll contains 
concluding remarks. 

II. THE MODEL 

We consider a large spin, S 3> 1, in a strong stationary 
magnetic field along the easy axis z. The spin is mod- 
ulated by a transverse magnetic field with frequency top 
close to the Larmor frequency loq. The Hamiltonian of 
the spin has the form 

H =u S z - \DS 2 Z - S x Acosuj F t (h=l) (1) 

This Hamiltonian well describes many single-molecule 
magnets, including Mni2 crystals; D characterizes the 
magnetic anisotropy and A is the modulation amplitude. 
It also describes a nuclear spin, with D characterizing the 
quadrupolar coupling energy to an electric field gradient 
in the crystal with an appropriate symmetry. 

We assume that the Zeeman energy levels in the ab- 
sence of modulation are almost equidistant. We also as- 
sume that the resonant modulation is not too strong. 
These conditions are met provided 

\u -Uf\,DS,A<z:uo. (2) 

For many SMMs the inequality DS <C ujq is fairly de- 
manding and requires strong static magnetic fields; for 
example D « 0.6 K for Fe 8 (where S = 10) 3 . On the 
other hand, for more isotropic SMMs the anisotropy is 
much smaller; for example, D « 0.04 K for Fei7 where 
S = 35/2, see Ref. HI (our definition of D differs by a 
factor of 2 from the definition used in the literature on 
SMMs) . The anisotropy is usually much weaker for large- 
st nuclei and the condition (|5J) is not restrictive. 

The quantum dynamics of an isolated spin with Hamil- 
tonian Eq. |T]) was considered earlier—. Here we are in- 
terested in the spin dynamics in the presence of dissi- 
pation. Different dissipation mechanisms are important 
for different systems. For single-molecule magnets, en- 
ergy dissipation is due primarily to transitions between 
spin energy levels accompanied by emission or absorp- 
tion of phonons. The transitions between both nearest 
and next nearest spin levels are importan t 18 ' 19 ' 23 . The 
corresponding interactions are 

H[ 1] = E k KS+S a + S Z S+) b k + h.c] (3) 
H\ 2) = E k Vt 2) (S\b k + h.c.) , s ± = s x ± l s y , 

where k enumerates phonon modes, b k is the annihilation 
operator for the k-th mode, and Vu 1 ^ and are the 
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coupling parameters responsible for transitions between 
nearest and next nearest Zeeman levels. The phonon 
Hamiltonian is 



H, 



ph 



(4) 



A similar interaction Hamiltonian describes the coupling 
of a nuclear spin to phonons, cf. Ref. and the early 
work 25 ' 26 . 

Along with the interaction ((3J we will consider the in- 
teraction that is linear in the spin operators, 



(5) 



Such interaction is allowed by time-reversal symmetry 
in the presence of a strong static magnetic field, with 

(3) 

the coupling constants V fe proportional to the off power 
of the field field. It can be thought of as arising from 
phonon-induced modulation of the spin g-factor. The 
interaction Eq. ([5]) is also important for impurity spins 
in magnetic crystals, in which case bk is the annihilation 
operator of a magno n 27 i 28 . 



A. Rotating Wave Approximation 

The dynamics of a periodically modulated spin can 
be conveniently described in the rotating wave approxi- 
mation. To do this we make a canonical transformation 
U(t) = exp(—iuj F S z t). The transformed Hamiltonian Hq 
then becomes H = U*H U - iWtl, 



H Q = -5uS z -\DS 2 z -\AS x , 

5ljJ = U)p — LU . 



(6) 



Here we disregarded fast-oscillating terms 
oc Aexp(±2iuj F t). 

We note that the Hamiltonian (j^J) has the form of a free 
energy of a magnetic moment in an easy axis ferromag- 
net, with S playing the role of the magnetization and 5lu 
and A giving the components of the effective magnetic 
field (in energy units) along the z and x axes, respec- 
tively. 

It is convenient to change to dimensionless variables 
and rewrite the Hamiltonian as Hq = S 2 D(g + p 2 /2), 
with 

1 
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s = S/S, 



^(s z +p) 2 - fs x , 



variables s are advantageous for describing large spins, 
since the commutators of their components are oc S , 
which simplifies a transition to the classical limit for 
S» 1. 



B. Quantum kinetic equation 

We will assume that the interaction with phonons 
(magnons) is weak. Then under standard conditions the 
equation of motion for the spin density matrix p is Marko- 
vian in slow time t, i.e., on a time scale that largely ex- 
ceeds LOp 1 and the typical correlation time of phonons 
(magnons) . We will switch to the interaction representa- 
tion with respect to the Hamiltonian uj f S z + -ff p h- Then 
to leading order in the spin to bath coupling the quantum 
kinetic equation can be written as 



s-V^kd-fWp-fWp-fw, 



(8) 



where A = dA/dr. 

The operators describe relaxation due to the in- 
teractions H> J> , with j = 1,2,3. They can be written 
schematically as 



Tp = r[(n+l) (L + Lp-2LpL + + pL+L) 
+n (LL+p - 2L+pL + pLL+)\ 



(9) 



Here we have taken into account that all transitions be- 
tween spin states with emission or absorption of phonons 
(magnons) involve almost the same energy transfer AE, 
with AE w u F for terms oc r (1) ,T (3) and AE w 2lj f 
for the term oc r*- 2 -* . In this sense, equation for spin 
relaxation §§§ resembles the quantum kinetic equation 
for a weakly nonlinear oscillator coupled to a bosonic 
bath 17 . Respectively, n is the Planck number of the 
emitted/ absorbed phonons, n = [exp(AE/kT — 1)] . 
Because of the same transferred energy, different transi- 
tions are characterized by the same rate constants, which 
for the interactions £T- ^ have the following form, in 
dimensionless time: 

r( 1 ) = ^- 1 5 2 ^Jy fe (1) | 2 5(c F -c fc ), 



=7TD- 1 S 2 J2 k \V k {2) S(2u F -u k ) 



H = 5u/SD f = A/2SD (7) The operators L for the interactions H, t 



(l)-(3) 



are 



The Hamiltonian g describes the dynamics of an isolated 
spin in "slow" dimensionless time r — SDt. It gives 
dimensionless quasienergies of a periodically modulated 
spin in the RWA. From Eq. {7J, the spin dynamics is de- 
termined by the two dimensionless parameters, p and /, 
which depend on the interrelation between the frequency 
detuning of the modulating field 5u>, the anisotropy pa- 
rameter DS, and the modulation amplitude A. The spin 



(10) 



(11) 



where s± = S±/S. 

It is important to note that, along with dissipation, 
coupling to phonons (magnons) leads to a polaronic ef- 
fect of renormalization of the spin energy. A standard 
analysis shows that renormalization due to H\ , to sec- 

(3) 

ond order in H\ comes to a change of the anisotropy 
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parameter D and the Larmor frequency. A similar change 
comes from nonresonant terms oc S + b^ k + H.c. In con- 
trast, renormalization from H^'^ 2 \ along with terms 
oc S z ,Sz, leads to terms of higher order in S z in the 
spin Hamiltonian, in particular to terms oc S z . The con- 
dition that they are small compared to the anisotropy 
energy DS Z imposes a constraint on the strength of the 

coupling ff} 1 ^ . Respectively, we will assume that the 
dimensionless decay rates r^ 1 ^ 2 ) are small, rW'( 2 ) <C 1. 
It is not necessary to impose a similar condition on the 
dimensionless rate I^ 3 ). Still we will be interested pri- 
marily in the spin dynamics in the underdamped regime, 
where r^ - ^ are small. 



III. CLASSICAL MOTION OF THE 
MODULATED SPIN 

The analysis of spin dynamics is significantly simplified 
in the classical, or mean-field limit. Classical equations of 
motion for the spin components can be obtained by mul- 
tiplying Eq. ([5]) by Sj (i = x,y,z), taking the trace, and 
decoupling Tr (s^s^p) — > s^s^. The decoupling should 
be done after the appropriate commutators are evaluated; 
for example, Tr ([s z ,g] p) — ► —ifs y . From Eqs. (O, 10), 
(fTTj) we obtain 

s = -s x d s g + (s) d , (s) d = r d (s 2 )s x (s x z) , 

T d { Sz ) = 2 (4r« s 2 + 2r< 2 )(i - si) + r< 3 >) , (12) 

where z is a unit vector along the z-axis, which is the 
direction of the strong dc magnetic field. 

We have assumed in Eq. (fT5|) that S » h. Note that, 
in dimensional units, S = |L|/?i, where L is the angular 
momentum, whereas in the classical temperature limit 
h = kTjfvjjp or kT/2hhjp depending on the scattering 
mechanism. Therefore the condition S 3> n imposes a 
limitation on temperature. 

Equation (fT2|) is reminiscent of the Landau-Lifshitz 
equation for magnetization of a ferromagnet. However, in 
contrast to the Landau-Lifshitz equation a retardation- 
free equation of motion for a classical spin could be ob- 
tained only in the rotating frame, that is, in slow time r. 
The term with d s g describes precession of a spin with en- 
ergy (quasienergy, in the present case) g. The term (s) d 
describes the effective friction force. It is determined by 
the instantaneous spin orientation, but its form is differ- 
ent from that of the friction force in the Landau-Lifshitz 
equation. 

We emphasize that Eq. (fT2|) is not phenomenological, 
it is derived for the microscopic model of coupling to the 
bath ©, ([5]). We now consider what would happen if 
we start from the Landau-Lifshitz equation and switch to 
the rotating frame using the RWA in the assumption that 
the resonant driving is comparatively weak, A<wo [cf. 
Eq. @]. In this case one should keep in the expression 
for the friction force only the leading term in the effective 



magnetic field, i.e., assume that H || z. The result would 
be Eq. (fT2")l with a dissipative term of the same form as 
the term oc but without dissipative terms that have 
the structure of the terms oc r^^r^ 2 ^. However, these 
latter terms play a major role for SMM a 18 ' 19 ' 20 ' 21 and for 
phonon scattering by nuclear spins. 

As mentioned in the Introduction, the dynamics of 
a single-domain magnetic nanoparticle in a circularly 
polarized field was studied using the Landau-Lifshitz- 
Gilbert equation in a series of papers^. It is clear from 
the above comment that the results of this analysis do not 
generally describe resonant behavior of SMMs. More- 
over, per iodic states in the rotating frame predicted in 
Ref. [jji do not arise in resonantly excited spin systems 
with relaxation oc , as shown below. 



A. Stationary states in the rotating frame for weak 
damping 

A classical spin is characterized by its azimuthal 
and polar angles, <fi and 9, with s z = cos 9, s x = 
sin 9 cos 4>, s y — sin 9 sin <fi. In canonically conjugate vari- 
ables (j>, s z equations of motion (fT2|) take the form 

4> = d Sz g, (13) 
Sz = -d^g - T d (s z )(l - s 2 ), 

where g as a function of s z , </> has the form g = — (s z + 
/i) 2 /2 - /(l - s 2 ) 1 / 2 cos0, cf. Eq. 0. We note that the 
dissipation term is present only in the equation for s z . 

In the absence of relaxation, precession of a spin with 
given g corresponds to moving along orbits on the (0, s z )- 
plane. The orbits can be either closed or open; in the 
latter case </> varies by 27r over a period, cf. Fig. [T] There 
are also stationary states where the spin orientation does 
not vary in time. Generally, relaxation breaks this struc- 
ture. If it is weak it makes some of the stationary states 
asymptotically stable or unstable and can also transform 
some of the orbits into stable or unstable limit cycles, 
which correspond to periodic oscillations of s z and <p in 
the rotating frame. The frequency of these oscillations 
is determined by the system nonlinearity and is not im- 
mediately related to a combination of the modulation 
frequency and the Larmor frequency, for example. 

Stationary states of Eq. (|13p , which is written in the ro- 
tating frame, correspond to the states of forced vibrations 
of the spin components s x , s y at frequency top in the lab- 
oratory frame. Periodic states in the rotating frame cor- 
respond, in the laboratory frame, to periodic vibrations 
of s z and to vibrations of s x , s y at combination frequen- 
cies equal to lof with added and subtracted multiples of 
the oscillation frequency in the rotating frame (which is 
small compared towj?). In what follows we keep this cor- 
respondence in mind, but the discussion refers entirely to 
the rotating frame. 

The analysis of stability of stationary states is based on 
linearizing the equations of motion near these states and 
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FIG. 1: (Color online). Phase portraits of the spin on (0,4>)- 
plane (s z = cos 61). The data refer to = T (2) = 0, T (3) = 
0.1, and / = 0.3. In panels (a)-(d) fx = -0.6,-0.2,0,0.2, 
respectively. 



looking at the corresponding eigenvalues Ai, A2 22 . In the 
absence of damping the stationary states are either hy- 
perbolic points (saddles) with real Ai^ or elliptic points 
(centers) with imaginary Ai,2- From Eq. (|13| . a fixed 
point is hyperbolic if A1A2 = T> < 0, where 



V = dlgdlg-(d 4> d Sz g) 2 



(14) 



(the derivatives are calculated at the stationary state). 
On the other hand, if V > the stationary state corre- 
sponds to an elliptic point, orbits g — const are circling 
around it. 

For weak damping, hyperbolic points remain hyper- 
bolic. On the other hand, a center becomes asymptot- 
ically stable (an attractor) or unstable (a repeller) for 
T < or T > 0, respectively. Here T = -d{T d (s z )(l ~ 
s 2 )]/<9s z , or in explicit form 



r = -At 



4r«(i - 2s 2 z ) - 4r( 2 >(i - s 2 z ) ~ r^l , (15) 



where s z is taken for the appropriate center; Ai + A2 = T . 
The sign of T determines stability of a stationary state 
also where dissipation is not small. 

The quasienergy g has symmetry properties that the 
change / — * — / can be accounted for by replacing 
<p —>■ cj) + 7r, s z — ► s z . This replacement preserves the 
form of equations of motion (fT3|) also in the presence of 
damping. Therefore in what follows we will concentrate 
on the range / > 0. On the other hand, the change 
fi — > — fi would not change g if we simultaneously replace 



4> — > 4>, s z — > — s z . In equations of motion one should ad- 
ditionally change r — * — r. Therefore, if for fi = fj,(°> < 
the system has an attractor located at a given (<^°) , s z °^), 

then for fi = —fiQ it has a repeller located at 4>^°\ —s z °\ 
This behavior is illustrated in Fig. [T] where panels (b) 
and (d) refer to opposite values of fi. 



B. Saddle-node bifurcations 

The function g(s) has a form of the free energy of a 
magnetic moment of an easy axis ferromagnet, as men- 
tioned earlier, with fi and / corresponding to the com- 
ponents of the magnetic field along and transverse to 
the easy axis, respectively. It is well known that g may 
have either two or four extreme points where dg/ds z = 
dg/d(f> = 0. The region where there are four extrema lies 
inside the Stoner-Wohlfarth astroid 29 |/| 2/3 + \fi\ 2/3 = 1 
on the plane of the dimensionless parameters fi and /, 
see Fig. [U[a) • The extrema of g outside the astroid are a 
minimum and a maximum, whereas inside the astroid g 
additionally has a saddle and another minimum or max- 
imum. All of them lie at 4> = or (f> = ir. 

In the presence of weak damping, the minima and max- 
ima of g become stable or unstable stationary states. We 
note that there are no reasons to expect that the sta- 
ble states lie at the minima of g, because g is not an 
energy but a quasienergy of the spin. The number of 
stable/unstable stationary states changes on the saddle- 
node bifurcation curve on the (/, /x)-plane. The condition 
that two stationary states merged 2 - has the form 



V + Td 4> d Sz g = 0. 



(16) 



For weak damping a part of the curve given by this 
equation is close to the astroid. On the astroid s z = 
— sgn(/i)|//| 1 / 3 . Then from Eq. (fT5)) for the merging sad- 
dle and node 



= -4sgn(/x) v /l-|/|2/3 

x (4r( 1 )(i-2|/| 2 / 3 ) + 4r( 2 )|/| 2 / 3 + r( 3 )) .(17) 



T 



If damping cx is weak, the node is stable for fi > 
and unstable for /i < 0. On the other hand, if is 
large, the stability depends on the value of /. 

The most significant difference between the saddle- 
node bifurcation curve and the Stoner-Wohlfarth astroid 
is that the bifurcation curve consists of two curvilinear 
triangles, that is, the astroid is "split" , see Fig. [2Kb) and 
Fig. [5] below. This is also the case for a modulated mag- 
netic nanoparticlei^. The triangles are obtained from 
Eqs. (fT5|) and (fT!))) . After some algebra we find that the 
"bases" of the bifurcation triangles are given by expres- 
sion 



/ s «±r d ( M )(i- M 2 ) 



2N1/2 



(18) 



to leading order in T^. This expression applies not too 
close to the vertices of the triangles. We note, however, 
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FIG. 2: (Color online). Saddle-node bifurcation lines. Panel 
(a): zero-damping limit, the lines have the form of the Stoner- 
Wolfarth astroid in the variables of reduced amplitude / and 
frequency detuning fi of the resonant field. Panel (b): nonzero 
damping, T (3) = 0.1, = T (2) = 0. In the dashed region 
the spin has two coexisting stable equilibria in the rotating 
frame. 



that Eq. (fT5|) gives the exact bifurcational value of /b for 
fi = and arbitrary 1^(0). 

The shape of the gap between the upper and lower 
curvilinear bifurcation triangles depends on the damping 
mechanism. In particular, the damping oc T^ 1 ' does not 
contribute to the gap for small (cf. Fig. O, whereas 
the damping oc does not contribute to the gap at 
small 1 — \fi\. The damping-induced change of the sides 
of the triangles compared to the astroid is quadratic in 
Td, far from the small-/ range. 

The positions of the small-/ vertices of the bifurcation 
triangles fc^c for small damping can be found from 
Eqs. (fT3]) and the condition that Eq. (fl"6|) has a degener- 
ate root, which gives 

f i c -±[l-V3(-T 2 d + TT d ) 1/2 ] , 
fc « ±(64/27) 1 / 4 T 3 / 4 (T d + (1/2)T) 1/2 (-r d + T) 1/4 
where and T are calculated for s z = 1. 

C. Periodic states and Hopf bifurcations 

An important property of the modulated classical spin 
is the possibility to have periodic states in the rotating 
frame. Such states result from Hopf bifurcations in which 
a stationary state transforms into a limit cycled. A Hopf 
bifurcation occurs if 



T = 0, 



V > 



in the stationary state. Besides the special case s z = 
discussed in Sec. IIV1 the corresponding stationary state 
is at s z = s z h, where 

i /4rW -4r< 2 ) -r( 3 )\ 1/2 

s *" = a 2\ 2TW - W) ) ' (19) 

a = ±i, r« > r< 2 ) + 

~ 4 



(the inequality on the damping parameters follows from 
the condition (s 2 )# < 1) 

The field Jh on the Hopf bifurcation lines as a func- 
tion of the reduced detuning /i is given by a particularly 
simple expression for weak damping. In this case, from 
second equation (fl3|) the phase <pn for the bifurcating 
stationary state is close to cither or tt with the addi- 
tional constraint d^gd^g > 0. Then from first equation 
(fi"3|) and Eq. ([19]) we find that Hopf bifurcation curves are 
straight lines, in the limit of vanishingly small damping, 



f H = ±[l-sl H ] 1/2 [1+m^] 



\f H \> [l-4 H f 2 or \v\>\s zH \ 



(20) 
(21) 



The structure of these lines is seen in Fig. \5\ below. They 
end on the saddle-node bifurcation curves and are tan- 
gent to these curves at the end points. A detailed analysis 
is presented in Sec. I VII 



HAMILTONIAN-LIKE MOTION AT EXACT 
RESONANCE 



IV. 



The spin dynamics (|12| displays an unusual and un- 
expected behavior where the modulation frequency top 
coincides with the Larmor frequency ujq, in which case 
/i = 0. This is a consequence of the symmetry of the 
quasienergy and the dissipation operator. In a certain 
range of dynamical variables 4>,s z , the spin behaves as 
if there were no dissipation, even though dissipation is 
present. This behavior is seen in the pattern of phase 
trajectories of the spin. An example of the pattern is 
shown in Fig. QTc) for the case T^ 2 ) = T^ 3 ) = 0, but 
the behavior is not limited to this case. As seen from 
Fig. [He), phase trajectories form closed loops, typical 
for Hamiltonian systems. 

For | /| lying inside the bifurcation triangles, the 
Hamiltonian-like dynamics occurs only in a part of the 
phase plane. This region of / corresponds to Td(0) < 
l/l < [1 + T^O)] 1 / 2 [the upper bound on |/| for fj, = 
can be easily obtained from Eqs. (fT5)) . (| 16[) ] . Here, 
the spin has four stationary states. For small |/z| two 
of them have small \s z \, s z w — /i/(l — /cos^>) where 
sin</> w —Td{0)/f. One of these states is a saddle 
point [(f> « — arcsin[Fd(0)//]] and the other is a focus 
[0«7r + arcsin[r d (O)//]]. 

For /j, = there occurs a global bifurcation, a homo- 
clinic saddle-saddle bifurcation (saddle loop—) where the 
separatrix coming out from the saddle goes back into it, 
forming a homoclinic orbit. Simultaneously, the focus in- 
side the loop becomes a center, T = for s z = 0. All 
trajectories inside the homoclinic orbit are closed loops. 
In contrast to the case of the vicinity of the double-zero 
eigenvalue bifurcation^ 2 ., the pattern persists throughout 
a broad region of /. 

We show how a Hamiltonian-like region in phase 
space emerges first for weak damping. For fi = the 
quasienergy g corresponds to the Hamiltonian of a spin 
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with anisotropy energy oc S 2 , which is in a transverse 
field oc /. Such spin in quantum mechanics has special 
symmetry, it can be mapped onto a particle in a symmet- 
ric potential2*i£. A part of the classical g — const orbits 
are closed loops on the (0, s z )-plane. They surround the 
center (s z = 0,4> = 7r). The orbits are symmetric with 
respect to the replacement 

s z -> ~s z , 4>^ 4>, (22) 

which leads to <p ~ * —4>i s z — * s z . 

Weak damping would normally cause drift of 
quasienergy. The drift velocity averaged over the period 
T p (g) of motion along the orbit is 

(.<?) = -r- 1 I " drd s ^gT d ( Sz )(l - s 2 z ). (23) 
Jo 

From the symmetry (|2"2")l and the relation Td(s z ) = 
Fd(~s z ), we have (g) = on a closed orbit for /i = 0. 
Therefore a closed orbit remains closed to first order in 
T e i- Of course, for open orbits, where <p is incremented by 
2n over a period, (g) ^ 0. These orbits become spirals in 
the presence of damping. 

Spirals and closed orbits should be separated by a sep- 
aratrix, which must be a closed orbit itself. Since the 
separatrix must start and end at the saddle point, we 
understand that at fi = for small Td there occurs a 
saddle-saddle homoclinic bifurcation. 

The topology discussed above persists as Td increases. 
The symmetry (|22[) is not broken by Td- Indeed, from 
equations of motion (|13[) . any orbit that crosses s z — 
twice per period for /i = has the property (j2"2"|) and 
therefore is closed. The closed orbits surround the center 
s z = 0, 4> = 7r — arcsin(rd(0)//) and fill out the whole 
interior of the separatrix loop. 

The Hamiltonian-like behavior is displayed also for fi = 
and / lying outside the bifurcation triangles. Here, 
the system has two stationary states, both with s z = 
but with different </>. From Eq. (fT5|) . for both of them 
T changes sign as fi goes through zero. Because there 
is no saddle point, for small \fi\ there is no separatrix, 
trajectories spiral toward or away from stationary states 
and possibly limit cycles. It follows from the arguments 
above that for /x = all trajectories become closed orbits. 
This is confirmed by numerical calculations for different 
relaxation mechanisms. 

It is convenient to analyze the overall dynamics of the 
spin system for /i / separately for the cases where 
the system does or does not have stable periodic states 
in the rotating frame. In turn, this is determined by 
the interrelation between the damping parameters, cf. 
Eq. (jTSJ) . Such analysis is carried out in Sees. fVl and IVII 

V. SPIN DYNAMICS IN THE ABSENCE OF 
LIMIT CYCLES 

We start with the case where the system does not have 
limit cycles. It corresponds to the situation where the 



damping parameter is comparatively small and the 
interrelation between the damping parameters (|19[) docs 
not hold. To simplify the analysis we set = T^ 2 ) = 0, 
i.e., we assume that the coupling to the bath is linear 
in the spin operators and is described by the interaction 
Hamiltonian H> ' . The qualitative results of this Section 
apply also for nonzero r^,r^ as long as +4T ( - 2 "> > 
4rW. The bifurcation dia gram for this case is shown in 

Fig. m 

From the form of the function T, Eq. (fTS]) , it follows 
that the damping oc T^ 3 ) transforms centers of conser- 
vative motion with s z > into unstable foci (repellers), 
whereas centers with s z < are transformed into stable 
foci (attractors). Therefore for /j, < the spin has one 
stable state. It also has one stable state in the unshaded 
region of the half-plane fi > (outside the bifurcation 
triangles in Fig. \2§- Inside the shaded regions within the 
triangles the spin has two coexisting stable states. 

Examples of the phase portrait are shown in Fig. [1] 
As expected, for weak damping the system has a stable 
and an unstable focus outside the bifurcation triangles, 
Fig. [Ha). I n t ne shaded region inside the triangle it has 
two stable foci, an unstable focus, and a saddle point, 
Fig. Did). In the unshaded region inside the triangle there 
is one stable and two unstable foci, Fig. [T](b) (the values 
of \x in panels (b) and (d) differ just by sign). 

A. Hysteresis of spin response in the absence of 
limit cycles 

The presence of two coexisting stable states leads to 
hysteresis of the spin response to the external field. Such 
hysteresis with varying dimensionless parameter /i, which 
is proportional to the detuning of the field frequency, is 
shown in Fig. [3J For large negative \x the system has one 
stable state with negative s z , cf. Fig.QJa). As /u increases 
the system stays on the corresponding branch (the lowest 
solid line in Fig. [3]) until the stable state merges with the 
saddle point (the saddle-node bifurcation). This happens 
for /i > as /i reaches the bifurcation triangle. As /i 
further increases the system switches to the branch with 
larger s z and then moves along this branch. If fi decreases 
starting with large values where the system has only one 
stable state, the switching to the second branch occurs 
for /i = 0. 

The hysteresis pattern in Fig. [3] differs from the stan- 
dard iS-shape characteristic. This is the case for any / 
lying between the minimum and maximum of the bi- 
furcation triangle for fj, — 0, i.e., for 2r^ 3 ^ < |/| < 

1/2 

(l + 4r' 3 ) 2 ) . It is a consequence of the fact that the 
bifurcation at /j, = is not a saddle-node bifurcation, 
whereas a most frequently considered S'-shapc hysteresis 
curve arises if both bifurcations are of the saddle-node 
type. In our case, for jj, = the branch which is sta- 
ble in the range of large positive fi (the upper stable 
branch in Fig. [3]) becomes unstable as a result of the 
motion becoming Hamiltonian-like. The value of s z on 
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this branch for [i = is s z = 0, it coincides with the 
value of s z at the saddle (but the values of s x are differ- 
ent). Therefore when s z is plotted as a function of /i the 
branch, which is stable for large positive fi crosses with 
the branch that corresponds to the saddle point. For neg- 
ative fi the branch, which is stable for large positive /i, 
becomes unstable, cf. Fig. [TJ As fj, decreases and reaches 
the bifurcation triangle for /i < 0, the saddle merges with 
an unstable equilibrium as seen in Fig. [3] 



1 

\ 

\ \ 

V. 









FIG. 3: (Color online). Hysteresis of spin response in the 
absence of periodic states in the rotating frame. The data 
refer to = T (2) = 0, T (3) = 0.1, and / = 0.3. The 

solid and dashed lines show, respectively, stable and unstable 
stationary states, the dotted line shows the saddle point. 

The spin components display hysteresis also if the 
shaded area of the bifurcation triangle in Fig. [H^b) is 
crossed in a different way, for example, by varying /. If 
the crossing occurs so that the line fi = is not crossed, 
the hysteresis curves have a standard S shape. We note 
that hysteresis of s x , s y corresponds to hysteresis of am- 
plitude and phase of forced vibrations of the spin. 



B. Interbranch switching without hysteresis 

The occurrence of Hamiltonian dynamics for fj, = 
leads to an interesting and unusual behavior of the sys- 
tem even outside the bifurcation triangles, i.e. in the 
region where the system has only one stable state. In 
the small damping limit and for |/| > 1 and \fj,\ <C 1 
the stationary states are at <j> = and <j> = tt, with 
s z = /V(/cos0 — 1)- The stable state is the one with 
s z < 0, whereas the one with s z > is unstable. As 
fi goes through zero the states with <j) — and (f> = tt 
interchange stability. This means that s x ~ cos 4> jumps 
between —1 and 1. Such switching is seen in Fig. 0] 



VI. SPIN DYNAMICS IN THE PRESENCE OF 
LIMIT CYCLES 

The classical dynamics of the spin changes significantly 
if the spin has stable periodic states in the rotating frame. 



FIG. 4: (Color online). Frequency dependence of the trans- 
verse spin component for field amplitudes / where the system 
has one stable state. The solid and dashed lines show the 
stable and unstable values of s x in the rotating frame. The 
data refer to = T (2) = 0, T (3) = 0.1, and / = 1.1. As the 
scaled frequency detuning /x goes through fj, = the value of 
s x changes to almost opposite in sign. 



This occurs where condition (|19p on the damping pa- 
rameters is met. The features of the dynamics can be 
understood by setting = = 0, > 0, i.e., by 
assuming that damping is due primarily to coupling to a 
bath H\ ; , which is quadratic in spin components, with 
elementary scattering processes corresponding to transi- 
tions between neighboring Zeeman levels. This model is 
of substantial interest for single-molecule magnet o 18 i 21 . 




(b) 



(i) 



(iii) 



(iv) / 


\ (ii) 
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FIG. 5: (Color online), (a) Saddle-node bifurcation lines for 
= 0.05, T (2) = T (3) = 0. (b) Saddle-node (solid lines) and 
Hopf bifurcation (dotted lines) in the limit of small damping 
oc r' 1 - 1 . Not too close to the astroid (see Sec. I VI B| for weak 
damping the system has the following states: (i) a stable and 
an unstable focus; (ii) two unstable foci and a stable limit 
cycle; (iii) a stable and an unstable focus and a stable and an 
unstable limit cycle; (iv) two stable foci and an unstable limit 
cycle. 

The saddle-node bifurcation curves for weak damping 
oc rW are shown in Fig. [5j Inside the curvilinear trian- 
gles the spin has four stationary states, whereas outside 
the triangles it has two stationary states. In contrast 
to the case of damping oc shown in Fig. [2j in the 
present case the bases of the triangles touch at n = 0. 
From Eq. (fTTjl . one of the states emerging on the sides 
of the triangles is stable for > 0, |/| < 2~ 3 / 2 and is 
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unstable otherwise; note that the stability changes in the 
middle of the bifurcation curves. 

The occurrence of periodic oscillations of the spin is 
associated with Hopf bifurcations. In the present case, 
from Eq. (flU)) the Hopf bifurcational values of s z are 
s zH = ±1/V2. Therefore Eq. QQ) for the Hopf bifur- 
cation lines for weak damping takes a simple form 

f H = 2- 1 ' 2 ±^ f H i (0,2- 3 / 2 ); (24) 
/* = -2-V»± M) f„i (-2-^,0). 

These lines are shown in Fig. \5\b). For |/| ~ 1 and far 
from the end points of the bifurcation lines, the typical 
frequency of the emerging oscillations is ~ 1 in dimen- 
sionless units, or ~ DS/h in dimensional units. 

A. Phase portrait far from the astroid 

Evolution of the spin phase portrait with varying pa- 
rameters far away from the astroid, \/j,\ 3> 1, can be un- 
derstood by looking at what happens as the Hopf bifur- 
cation curves are crossed, for example by varying /. The 
result is determined by two characteristics. One is stabil- 
ity of the stationary state for / close to the bifurcational 
value fjj . The stability depends on the sign of T for small 
/ — Jh (note that T changes sign for / = fjj). The other 
characteristic is the sign of the quasienergy drift velocity 
(g) for / = fa and for g close to its bifurcational value 
gn at the stationary state. It is given by Eq. (|23p [note 
that, generally, (g) oc (g - g H ) 2 for / = f H \- A combi- 
nation of these characteristics tells on which side of the 
bifurcation point there emerges a limit cycle and whether 
this cycle is stable or unstable. 

We write the value of s z at the Hopf bifurcation 
point as s z h — a/y2, where a — ±1, cf. Eq. (fT!))) . 
The bifurcational value of the field (|24| is fjj = 
± (2 -1 / 2 + a/i) cos 4>h , where 4>h is the phase of the 
bifurcating stationary state. Linearizing Eq. (|15p in 
s z — s z h and using the explicit form of the determinant 
T> one can show that, for small / — in a stationary 
state sgn[T/(/ - f H )} = -sgn.[af H ]- Then 

sgnT = -{a sgn/ ff ) sgn(/ - f H ). (25) 

The analysis of the quasienergy drift velocity near a 
Hopf bifurcation point is done in Appendix It follows 
from Eqs. ([AT]) . ([X2|) that 

(g) = CaY^(g - g H ) 2 (j3\f H \ - Vt) , 
sga[(g)/(g - g H )) = a/? sgn (p\f H \ - yty , (26) 

where C > is a constant and j3 = sgn(/# cos (j>H ) = 
sgn(2" 1 / 2 + afjt) = ±1 [fi is related to f H by Eq. ([M]) ]. 

The sign of (g)/{g — gn) shows whether g approaches 
gn as a result of damping or moves away from g^. If 
sgn [(g) (g — gn)] < 0, the vicinity of the stationary state 



and the nascent limit cycle attracts phase trajectories. 
Therefore at a Hopf bifurcation a stable focus becomes 
unstable and a stable limit cycle emerges. On the other 
hand, if sgn [(g)(g — gn)\ > 0, at a Hopf bifurcation an 
unstable focus transforms into a stable one and an un- 
stable limit cycle emerges. 

Equation (|25p allows one to say on which side of 
i.e., for what sign of / — fn the stationary state is sta- 
ble, since for a stable state T < 0. Therefore together 
Eqs. (|25|) and (|26|) fully determine what happens as / 
crosses the bifurcational value. 

We are now in a position to describe which states ex- 
ist far from the astroid in different sectors (i)-(iv) in 
Fig. [3(b). For small |/| and large regions (i) in 
Fig. [5^b) , the system is close to a spin in thermal equi- 
librium, it has one stable and one unstable stationary 
state. We now start changing / staying on the side of 
large positive fj,. When / crosses one of the bifurcation 
curves fn — ± (2 -1 / 2 — //) , the system goes to one of the 
regions (ii) in Fig. [5jb) . On the both bifurcation curves 
a = fi = — 1. Therefore, from Eqs. (|2"5|) , (f2"6")h when one 
of these curves is crossed as |/| increases, there emerges a 
stable limit cycle, and the stable focus becomes unstable. 
As |/| further increases it crosses the bifurcation curves 
±(2 -1 / 2 + fi) and the system goes to one of the regions 
(hi) in Fig. 0(b) (we assume that the crossing occurs in 
the region \fn\ > 2 1 / 2 ). On these bifurcation curves 
a — (3 = 1. Therefore, from Eqs. ([23)) . (j2l)j) . when they 
are crossed with increasing |/| there emerges an unstable 
limit cycle and the unstable focus becomes stable. 

We now start from the range of large negative /u, and 
small |/|. As we increase |/| and cross the bifurcation 
curves fn = ±(a* + 2 -1 / 2 ) the system goes from region 
(i) to one of the regions (iv) in Fig.[5]Jb). From Eqs. (|2"5j) . 
(f2"6")l . in this case an unstable focus goes over into a stable 
focus and an unstable limit cycle emerges. Further cross- 
ing into one of the regions (hi) with increasing |/| leads 
to a transformation of a stable focus into an unstable fo- 
cus and an onset of a stable limit cycle. These arguments 
were used to establish the nomenclature of states in re- 
gions (i)-(iv) in Fig. [5{b). They agree with the results of 
direct numerical calculations. 



B. Other bifurcations of limit cycles 

1. Merging of stable and unstable limit cycles 

The number of periodic states in the rotating frame 
may change not only through Hopf bifurcations, but also 
through other bifurcations, where the radius of the bi- 
furcating limit cycle does not go to zero. The simplest 
is a bifurcation where a stable limit cycle merges with 
an unstable limit cycle (saddle-node bifurcation of limit 
cycles). The onset of such bifurcations is clear already 
from Eq. (j2l)|) . Indeed, at a Hopf bifurcation point the 
equation for the period-averaged quasienergy has a form 
(9) = c(g-g H ) 2 + - ■ ■ with c cx 0[f H [~V2. For \f H \ = V2 
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on the bifurcation curves ([24]) with (3 = 1 [the top and 
bottom dotted lines in Fig. Efb)] the coefficient c = 0. 
This is a generalized Hopf bifurcation^, see Fig. |6] 

At the generalized Hopf bifurcation, in phase space 
(4>,s z ) a stationary state merges simultaneously with a 
stable and an unstable limit cycle. In parameter space 
(//,/), the Hopf bifurcation curve coalesces with the 
curve where stable and unstable limit cycles are merg- 
ing, and the latter curve ends. The bifurcation curves are 
tangent, the distance between them scales as a square of 
the distance to the end point /3\fn\ = V2 if the latter 
distance is small. This is seen in Fig. [5] In the compara- 
tively narrow region between the Hopf bifurcation curve 
and the corresponding limit-cycle merging curve the sys- 
tem has three limit cycles. One of these cycles disappears 
on the Hopf bifurcation curve, so that in regions (iii) in 
Fig. [3h) there are two limit cycles and deeper in regions 
(ii) and (iv) there is one limit cycle. On its opposite 
end, the curve of merging limit cycles coalesces with the 
saddle-loop bifurcation curve. 




FIG. 6: (Color online). Bifurcation diagram in the limit 
r' 1 ' — > 0. The diagram is symmetric with respect to /i = 
and / = axes, and therefore only the quadrant / > 0, fi < 
is shown. Saddle-node, Hopf, and saddle-loop bifurcation 
curves are shown by the solid, dotted, and long-dash lines, 
respectively, whereas the short-dash line shows the curve on 
which stable and unstable limit cycles merge. 



2. Saddle loops 

Spin dynamics for damping oc rw is characterized also 
by global bifurcations of the type of saddle loops. This 
is clear already from the analysis of the end points of the 
Hopf bifurcation curves. These points lie on the curves of 
saddle-node bifurcations. The corresponding equilibrium 
point has double-zero eigenvalue, and the behavior of the 
system near this point is well-known^. The Hopf bifur- 
cation curve is tangent to the saddle-node bifurcation 
curve at the end point. In addition, there is a saddle- 
loop bifurcation curve coming out of the same end point 



and also tangent to the saddle-node bifurcation curve at 
this point. At a saddle- loop bifurcation the system has a 
homoclinic trajectory that starts and ends at the saddle 
point. 

The structure of vicinities of the end points of the Hopf 
bifurcation curves is shown in Figs. [5] and [7] for the curves 
ending on the sides and the bases of the saddle-node 
bifurcation triangles, respectively. Note that the Hopf 
bifurcation curves that crossed at / = in the limit 
— > are separated for finite r^. They end on the 
saddle-node bifurcation curves. 




FIG. 7: (Color online). Bifurcation diagram near the end 
point of the Hopf bifurcation line which in the limit F' 1 - 1 — > 
has the form fn = —fi — 2 -1 ^ 2 . For nonzero r' 1 ' this 
bifurcation line ends on the saddle-node bifurcation line (|18|) . 
The plot refers to F^ 1 ' = 0.0125. The inset shows a close 
vicinity of the end point. Hopf, saddle-node, and saddle loop 
bifurcation curves are shown by dotted, solid, and long-dashed 
lines, respectively. Other Hopf bifurcation curves that go to 
fn = for r' 1 ^ — > display a similar behavior near their end 
points. 

We have found numerically a fairly complicated pat- 
tern of saddle- loop bifurcation curves. Full analysis of 
this pattern is beyond the scope of this paper. 

C. Hysteresis of spin response in the presence of 
limit cycles 

Coexistence of stable stationary states and stable limit 
cycles in the rotating frame leads to hysteresis of the re- 
sponse of a spin when the modulating field parameters 
are slowly varied. Examples of such hysteresis with vary- 
ing scaled frequency detuning fj, and the characteristic 
phase portraits are shown in Fig [5] 

The hysteretic behavior is unusual. This is a conse- 
quence of the feature of the spin dynamics for fi = 
where either all phase trajectories are closed loops (for / 
outside the curvilinear saddle-node bifurcation triangles 
in Fig. [5]) or all trajectories in a part of the phase plane 
are closed loops (for / inside the triangles in Fig. [5} . As 
a result two or more states (stationary or periodic) si- 
multaneously loose or acquire stability as fi goes through 
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(c) (d) 




FIG. 8: (Color online). Panels (a) and (b): hysteresis of the 
spin dynamics with varying scaled detuning of the modulating 
field frequency fj,. In (a) / = 0.4, so that fi goes through the 
curvilinear bifurcation triangle in Fig. [5] In (b) / = 1.2, it 
lies above the triangles. Bold solid lines, long dashed lines, 
and dotted line show stable and unstable equilibria and the 
saddle stationary state, respectively. Pairs of thin solid lines 
and short dashed lines show the boundaries (with respect to 
s z ) of stable and unstable limit cycles. Panels (c) and (d): 
phase portraits for /i = 0.2. In (c) and (d) / = 0.4 and 1.2, 
respectively. The arrows show the direction of motion along 
the trajectories. The data refer to r' 1 ' = 0.05. 

0. This leads to an ambiguity of switching, a "Buridan's 
ass" type situation. Where a stable branch looses sta- 
bility for fj, = 0, the system has more than one stable 
state to switch to. Also, in contrast to the situation of 
Sec. |V] where the system had no limit cycles, hysteresis 
emerges whether the varying field parameter does or does 
not cross the saddle- node bifurcation lines. 

Figures [5fa) and (b) show the behavior of the system 
with varying /i for / inside and outside the saddle-node 
bifurcation triangles, respectively. It should be noted 
that we chose / in Fig.|S^a) so that the saddle-loop bifur- 
cation line is not encountered, which provides an insight 
into the most basic features of the hysteresis. In addition, 
in Fig. EJb) we do not show an extremely narrow region 
near Hopf bifurcation lines /x « ±(/ — 2~ 1 / 2 ) where the 
system has small-radii stable and unstable cycles which 
merge on the short-dash bifurcation line in Fig. [S] 

In Fig. IDJa), for large negative \i the system has one 
stable state (with negative s z ). As \x increases this state 
disappears via a saddle-node bifurcation and the system 
switches to a stable limit cycle. For chosen / = 0.4 this 
happens for /i sa 0.33. With further increase of \x the 
limit cycle shrinks and ultimately disappears via a Hopf 



bifurcation, and then the stationary state inside the cycle 
becomes stable. 

On the other hand, if we start in Fig. [5fa) from large 
positive [i and decrease /x, the stable stationary state via 
a supercritical Hopf bifurcation becomes a stable limit 
cycle. The cycle looses stability at /x = 0, and as /i 
becomes negative the system can switch either to the 
stable stationary state inside the cycle (with s z — > +0 
for /i — > —0) or to a stable stationary state outside the 
cycle with negative s z . The stable state with s z — * +0 
for [i — > —0, ultimately looses stability with decreasing 
/i via a Hopf bifurcation (at /x » — / — 2 -1 / 2 , for small 
damping, cf. Fig. [5]). If the system is in this state, it 
switches to the stable equilibrium with negative s z . 

A typical phase portrait for / = 0.4, < \i < 0.33 is 
shown in Fig. [S£c) . It gives an insight into the behavior 
described above. The system has a stable limit cycle with 
an unstable focus inside and with stable and unstable 
equilibria and a saddle point outside the limit cycle. For 
fj, = the system has a homoclinic saddle connection, 
and all trajectories inside the homoclinic trajectory are 
closed loops, cf. Fig. QJc) 

In Fig. [5Jb), for large negative /x the system also has 
one stable state (with negative s z ). As /x increases this 
state looses stability via a Hopf bifurcation (at /x ~ 
— / + 2~ 1//2 , for small damping). The emerging state 
of stable oscillations looses stability for fi = 0. For larger 
/x the system switches either to the stationary state in- 
side the limit cycle (with s z — > +0 for /x — > +0) or to 
another stable periodic state. The coexistence of stable 
and unstable limit cycles with stationary states inside of 
them is seen in Fig. [8jd) . 

As /x becomes positive and further increases, the stable 
stationary state inside the unstable cycle looses stability 
by merging with this cycle, and the system switches to 
the periodic state corresponding to the stable limit cycle 
in Fig. [5Jd). For still larger /x (/x w / + 2 -1 / 2 , for weak 
damping) this state experiences a Hopf bifurcation and 
becomes a stable stationary state. The behavior with /x 
decreasing from large positive values can be understood 
from Fig. [5] in a similar way. 



VII. CONCLUSIONS 

We have developed a microscopic theory of a reso- 
nantly modulated large spin in a strong static magnetic 
field and studied spin dynamics in the classical limit. 
We have taken into account relaxation processes impor- 
tant for large-spin systems of current interest. They 
correspond to transitions between neighboring and next- 
neighboring Zeeman levels with emission or absorption of 
excitations of a bosonic thermal bath. Classical spin dy- 
namics depends significantly on the interrelation between 
the rates of different relaxation processes. Generally it is 
not described by the Landau-Lifshitz equation for mag- 
netization in a ferromagnet, although one of the coupling 
mechanisms that we discuss leads to the Landau-Lifshitz 



12 



damping in the rotating frame. 

We found that the spin dynamics has special symmetry 
at exact resonance where the modulation frequency is 
equal to the Larmor frequency, ujf = u>q. This symmetry 
leads to a Hamiltonian-like behavior even in the presence 
of dissipation. In the rotating frame, phase trajectories 
of the spin form closed loops in a part of or on the whole 
phase plane. Therefore when ujp goes through lu several 
states change stability at a time. 

The simultaneous stability change leads to unusual ob- 
servable features. Where the system has only one stable 
state for a given parameter value, as ujf goes through ojq 
there occurs switching between different states that leads 
to an abrupt change of the magnetization. The behav- 
ior is even more complicated where several stable states 
coexist for ujf close but not equal to ujq. Here, where 
u; f — changes sign, the state into which the system 
will switch is essentially determined by fluctuations or 
by history (if ujp is changed comparatively fast). 

We found the conditions where the spin has more than 
one stable stationary state in the rotating frame. Such 
stable states correspond to oscillations of the transverse 
magnetization at the driving frequency in the laboratory 
frame. Multistability leads to magnetization hysteresis 
with varying parameters of the modulating field. If the 
fastest relaxation process is transitions between neigh- 
boring states due to coupling quadratic in spin opera- 
tors, the resonantly modulated spin can have periodic 
nonsinusoidal states in the rotating frame with frequency 
oc DS/H, where D is the anisotropy energy. In the labora- 
tory frame, they correspond to oscillations of the trans- 
verse magnetization at combinations of this frequency 
(and its overtones) and the Larmor frequency. 

Quantum fluctuations of the spin lead to phase diffu- 
sion of the classical periodic states in the rotating frame. 
As a result, classical oscillations decay. The intensity 
of quantum fluctuations and the related decay rate de- 
pend on the value of S~ x . We have found^ that the 
oscillations decay comparatively fast even for S = 10. 
Therefore they are transient. Still the classically stable 
vibrations lead to pronounced features of the full quan- 
tum spin dynamics. 

The present analysis can be immediately extended to a 
more general form of the spin anisotropy energy, in par- 
ticular to the case where along with DS Z this energy has 
a term E(S 2 — S 2 ), which is important for some types of 
single-molecule magnets^. In the RWA, the correspond- 
ing term renormalizes D and ujq. The analysis applies 
also to decay due to two-phonon or two-magnon coupling, 
which often plays an important role in spin dynamics and 
leads to energy relaxation via inelastic scattering of bath 
excitations by the spin. Another important generaliza- 
tion is that the results are not limited to linearly polar- 
ized radiation. It is easy to show that they apply for an 
arbitrary polarization as long as the radiation is close to 
resonance. 

In conclusion, starting from a microscopic model, we 
have shown that the classical dynamics of a resonantly 



modulated large spin in a strong magnetic field dis- 
plays several characteristic features. They include abrupt 
switching between magnetization branches with varying 
parameters of the modulating field even where there is 
no hysteresis, as well as the occurrence of hysteresis and 
an unusual pattern of hysteretic inter-branch switching. 
These features are related to the Hamiltonian-like be- 
havior of the dissipative spin for modulation frequency 
equal to the Larmor frequency in the neglect of spin 
anisotropy. Along with forced vibrations at the mod- 
ulation frequency, the transverse spin components can 
display transient vibrations at a combination of the mod- 
ulation frequency and a slower frequency oc DS/h and its 
overtones. They emerge if the fastest relaxation mecha- 
nism corresponds to transitions between neighboring Zee- 
man levels with the energy of coupling to a thermal bath 
quadratic in the spin operators. 
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research was supported in part by the NSF through grant 
PHY-0555346 and by the Institute for Quantum Sciences 
at MSU. 



APPENDIX A: ENERGY CHANGE NEAR A 
HOPF BIFURCATION 

In this Appendix we outline the calculation of the re- 
laxation of quasienergy g near a Hopf bifurcation point. 
For concreteness we assume that T^ 2 ' = = and 
the only nonzero damping parameter is T^ 1 ). For small 
damping a stationary state that experiences a bifur- 
cation has phase 4>h close to either or 7r, whereas 
s z h ~ ±2~ 1//2 . The dynamics is characterized by two 
parameters, a = sgns z # and j3 = sgn[/# cos 4>h\- The 
bifurcational value of the field for T^ 1 ^ — > is fn = 

(2- 1 /2 +a/i ) CQS ff [ c f. Eq . (gj)]. 

At the bifurcating stationary state the quasienergy is 
gn = g {4>h, s z h); it is easy to see that this is a local 
minimum of g(0, s z ) for f3 > or a maximum for [3 < 0. 
Respectively, on phase plane (</>, s z ) the constant-g tra- 
jectories close to the bifurcating stationary state rotate 
about this state clockwise for j3 > and counterclock- 
wise for (3 < 0. The angular frequency of this rotation is 

2Tr/T p (g H ) = V 1 / 2 , where V is given by Eq. (fLf) . 

We now consider dissipation-induced drift over 
quasienergy (g). It is given by Eq. (|23j) . Noting that 
d s ,g = 4> an< i using the Stokes theorem we can rewrite 
this equation as 

(9) =Pr- x (g) J d<f>ds z T, (Al) 

where the integral is taken over the interior of the 
constant-g orbit on the (<fi,s z ) plane and T = T(s z ) is 
given by Eq. (|15p . At a Hopf bifurcation point T = 0. 
Therefore T(s z ) in Eq. (|A1[) must be expanded in Ss z — 
s z - s zH - 
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It is convenient to calculate integral (|A1[) by changing 
to integration over action-angle variables (I,ip), which 
are canonically conjugate to (s z , 4>), with g being the ef- 
fective Hamiltonian. The angle ip gives the phase of os- 
cillations with given quasienergy g. The action variable 
I = (2tt)~ 1 § s z d<j) is related to g by the standard expres- 
sion dl/dg — T p (g) /2n w I?" 1 / 2 ; we note that / becomes 
negative away from the stationary state for (3 < 0. 

In evaluating expression (|A1[) it is further convenient 
to start with integration over ip. It goes from to 2tt 
and corresponds to period averaging for a given I (X Sg = 
g — qh (integration over / corresponds to integration over 

Sg). 

If vibrations about (<Ph,s z h) were harmonic, the 
lowest-order term in 8s z that would not average to zero 
on integration over %j) would be (d 2 T ' J ' ds 2 z ){5s z ) 2 /2 oc \Sg\ 



(the derivative of T is calculated at the bifurcating sta- 
tionary state). However, it is easy to see that the integral 
over ip of the linear in 5s z term in T is also ~ Sg. It can 
be calculated from equation of motion (f> = d Sz g by ex- 
panding the right-hand side to second order in 6s z ,S<p 

and noting that = 0, where the over line means averag- 
ing over if). This gives, after some algebra, 

T= 64T^a(5g)(2 3 / 2 p\f H \-iy 2 

x(/3\f H \~2 1 / 2 ). (A2) 

This expression combined with Eq. (|A1 1) shows how the 
energy relaxation rate depends on the field fjj. It is used 
in Section IVll to establish the full bifurcation diagram. 
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